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Abstract 

The overdamped Brownian motion of a self-propelled particle which is driven by a projected 
internal force is studied by solving the Langevin equation analytically. The "active" particle 
under study is restricted to move along a linear channel. The direction of its internal force is 
orientationally diffusing on a unit circle in a plane perpendicular to the substrate. An additional 
time-dependent torque is acting on the internal force orientation. The model is relevant for active 
particles like catalytically driven Janus particles and bacteria moving on a substrate. Analytical 
results for the first four time-dependent displacement moments are presented and analysed for 
several special situations. For vanishing torque, there is a significant dynamical non-Gaussian 
behaviour at finite times t as signalled by a non-vanishing normalized kurtosis in the particle 
displacement which approaches zero for long time with a 1/t long-time tail. 
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I. INTRODUCTION 



The Brownian motion of self-propelled ( "active" ) particles bears much richer physics 
than the traditional diffusive dynamics of passive particles. Active particles can be mod- 
elled by moving under the action of an internal force sometimes combined with an internal 
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ace [111 ■ ■ 



or external torque. Realizations in nature are certain bacteria 
zoa [s], y, [l^ which swim in circles when confined to a sur 
is possible to prepare catalytically driven Janus particles 

cles \lm which perform self-propelled Brownian motion. For a recent investigation including 

n n 

confinement see pJ7|. On the macroscopic scale vibrated polar granular rods [18] on a planar 
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. S, y, 3| and spermato- 
n the colloidal world, it 



15j or biometric parti- 



substrate and even the trajectories of completely blinded and ear-plugged pedestrians pi 
can be considered as rough realizations of self-driven Brownian particles. If the particle is 
embedded in a liquid (a "swimmer"), as characteristic for colloids, the direction of its driv- 



ing force is fiuctuating, in general, according to orientational Brownian motion [20|,121|,122|. 
This gives rise to a non-ballistic translational motion of the particles which is coupled to the 
fiuctuating orientational degree of freedom. 

In most cases the direction of the self-propelling force is within the plane of motion. For 
colloidal particles, however, it is possible to confine the particle on a substrate by using, e.g.. 



strong gravity such that the particles are still freely rotating |15l . |23| . |24| though they are 
confined in a planar monolayer. In this situation the component of the self-propelling force 
which is normal to the surface is compensated by the substrate, i.e., only the projection of 
the self-propelling force onto the plane is driving the particle. Therefore the translational 
motion is coupled to the (Brownian) orientational motion 25 1. 



In this paper we consider a one-dimensional model [26| for the Brownian dynamics of a 
self-propelled particle on a substrate. The particle is self-propelled along its orientational 
axis, which itself is subjected to Brownian orientational diffusion. The particle is confined to 
a channel, however, such that only the projected force in channel direction is ac ting to drive 
the particle. The present study is more general than earlier work in reference 251]: first of 
all, the present calculation resolves the Cartesian components of the isotropic model on an 
unconfined plane. Second, an arbitrary time- dependence of the external torque is included 



here while this torque was constant in 



Finally, we calculate time-dependent moments 



of the particle displacement up to fourth order as compared to results up to second order 
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in reference 



251]. The results are discussed for several special cases. In general, long-time 



self- diffusion is found. Non-Gaussian behaviour is found for intermediate times as signalled 
in the corresponding fourth cumulant. The normalized kurtosis is positive for small times, 
then changes sign and approaches zero from below at long times with a 1/t long-time tail. 
This can be compared to recent investigations for an undriven Brownian ellipsoid [27^. In 
the latter case, the kurtosis was found to be positive approaching zero from above for long 
times with the same 1/t long-time tail. 

This work represents a first step towards a many-body situation of interacting self- 
propelled particles. These are also realizable in experiments (see, e.g., 12i_l^, Q). The 
suitable theoretical framework is the many-body Smoluchowski equation 2J, from which 
one can derive a coupled hierarchy of equations for the set of many-body distribution func- 



tions similar in s 
hierarchy 
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3irit to the traditional BBGKY (Bogolyubov-Born-Green-Kirkwood-Yvon) 
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3l| for Liouville dynamics, see also Felderhof [32] for a discussion in the 



context of Brownian motion. Therefore we think that this paper is particularly appropriate 
for this issue dedicated to the 100th anniversary of Prof. N. N. Bogolyubov. 

This paper is organized as follows: In section [Tll we propose and motivate the model. The 
first four displacement moments are calculated analytically for the torque-free case in section 
llllt while section IIVI contains the results for a general time-dependent torque. Finally, in 
section |Vl we conclude and give an outlook on possible future activities. 



II. THE MODEL 



The model system under study consists of a self-propelled colloidal sphere of radius 
i?, which is confined to an infinite linear channel in the x-direction, where it undergoes 
completely overdamped Brownian motion (for a sketch see figured]). Whereas the motion 
of the center-of-mass position x is constrained to one dimension, the orientation vector 
u = (cos</), sin0, 0) is constrained to rotate in the xy-plane. The self-propulsion of the 
particle is modelled by a constant effective force along the particle orientation F = Fqu and a 
generally time-dependent effective torque in the z-direction M = Me^. Because the particle 
is confined, only the projected force F ■ = FQC0S(j)ex drives the particle systematically 
along the channel. Based on these considerations, the translational and orientational motion 
is modelled by a Langevin equation for the center-of-mass position x and the orientation 
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FIG. 1: Sketch of the model system: A spherical colloidal particle (dark grey) is confined to a linear 
channel (light grey) along the x-direction. The self-propulsion is modelled by a constant effective 
force Fq along the particle orientation u. The latter is constrained to rotate in the xy-plane. Only 
the projected force FqcoscJ) drives the particle along the channel. A systematic, time-dependent 
torque M{t) = M{t)ez is also indicated. 



vector u: 



CxX 

— = (3D[FoCos(j) + f{t)] , 

— = f3D,[M{t)+g{t)]e, x u, 



(1) 
(2) 



where f{t) is a zero-mean, Gaussian white noise random force, which is characterized by 
(/(t)) = and {f{t)f{t')) = 25{t - t')/{/3^D), where angular brackets denote a noise aver- 
age. Correspondingly, g{t) is a Gaussian white noise random torque with {g{t)) = and 
{g{t)g(t')) = 26(t — t')/ {P^Dj.). Here, = k^T denotes the thermal energy. D and are 
the translational and rotational short-time diffusion constants, respectively. For a sphere of 
radius R in the three-dimensional bulk the two quantities fulfill the relationship 

D _ 4i?2 



(3) 



Due to the constraint on the orientational motion, the vector equation ([2]) reduces to a 
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Langevin equation for the orientational angle 0, which is given by 

^ = /3A mt) + g{t)] . (4) 

If the initial time to is set to be zero, the solutions of the Langevin equations ([1]) and (j4]) 
are given by 

0(t) = (3D, [ [M{t') + g{t')] dt' + <f)o (5) 
Jo 

and 

x{t)=(3D Fo f cos0(t')dt'+ / f{t')dt' +xo (6) 

L -'0 Jo \ 

with 00 = 'Pito) and Xq = x(to). 

The translation-rotation-coupling between these two equations, which is due to the cosine 
in equations ([1]) and (Q, leads to nontrivial results for the mean position {x — Xq) and the 
mean square displacement {{x — xq)^) of the particle position, as is shown in the following 
sections. Furthermore, the presence of the coupling term leads to non-Gaussian behaviour, 
which is reflected in a non-zero kurtosis. The latter is obtained by calculating the fourth 
moment of the particle displacement distribution further down. 

We start our analysis in section [TTll by studying the special case of a vanishing systematic 
torque M = 0. The more complex situations of a constant torque M(t) = M and a generally 
time-dependent torque Mit) are considered in section llVi 



III. RESULTS FOR A VANISHING TORQUE 

In this section, the simplest case with a vanishing torque is covered. Solving equation (jlj) 
for M(t) = and averaging gives 

m)) = 00 (7) 

and for the second moment 

{m)~M')=2D,t. (8) 



As 0(t) is a linear combination of Gaussian variables g{t'), according to Wick's theorem 20|, 
(f){t) is Gaussian as well. Thus the probability distribution of proves to be 
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FIG. 2: Mean position of a spherical particle without external torque for I3RFq = 10 and different 
values of . 



Now the mean position of the particle can be calculated. From 

/oo 
cos(0)P(0,t)d0 = e"-^^* COS 00 
-oo 

follows 

- xo) = ^l3FoR^ cos(0o) [l - e"^^*] , 
where we made use of equation ([3]). Thus for short times one obtains 

- xo) = ^/?Foi?2 cos(0o) At + O {e) 
and for t ^ Z^jT^ the 0o-dependent mean position converges towards 



lim (x(t) — xq) = -PFqR^ cos{(f)o)- 

t^oo 3 



(10) 



(11) 



(12) 



(13) 



The trajectory of the mean position (x(t)) is shown in figure [2] where the time t is given in 
units of -DjT^ while the length x is scaled by the particle radius R. 

To calculate the mean square displacement, the following integrals have to be solved: 



((x(t)-Xo)^> = P'D 



/?2 
-^0 



t pt 

dti / dt2(cos 0(ti) cos 0(t2)) 
Jo 

t ft ft ft 

+2Fo / dti / dt2(cos0(ti)/(t2))+ / dti / dt2(/(ti)/(t2)) .(14) 
Jo Jo Jo 



The third summand can be calculated easily and equals 2t/ (P'^D). As (cos 0(t)) only depends 
on the random torque g(t), (cos0(t)) and f(t) are statistically independent. Therefore 
the second summand vanishes. To calculate the first summand in equation (1141) . the time 
correlation function is used. With (pi = (j){ti) and 02 = ^(^2) the required average can be 
written as 

(cOS0iCOS02)t^>t2 = y #1 j d(p2COS(piCOS(p2G{(j)i,(i)2,tl-t2)P{<p2,t2)\ti>t2- (15) 

Here, G(0i, 02, — ^2) is the Green function, which is given by 

1 / f0i - 



^^^-^-'^-'^^ = V4.D.(t,-t2)^"" V^dJM^)) • ^^'^ 
This yields 

(cos 01 COS 02),^ = ie-^^(*^-*^) [1 + cos(20o)e-^^-*^] . (17) 

The expression for (cos0i cos02)j2>ti is obtained in exactly the same way by replacing ti 
and ^2 with each other. Now, the first summand in formula (fT4|) is calculated by simple 

integration and the mean square displacement can be written in the final form 

2 

e 

12 — - ^ 

(18) 

The long-time diffusion coefficient D\ is given by 



{{x{t)-x,f)=2Dt+[^l3FoR^^ 



A = ^hm i {{x{t) -x,Y) = D + ^ {f3FoR^yD,. (19) 

Figure [3] displays the results for the same cases that were examined in figure El The graph 
for 00 = TT coincides with the graph for 0o = 0. As can be seen in the logarithmic plots 
and from the expression f|T9l) . the initial angle 0o is not relevant for times much longer than 

In the following the non-Gaussian behaviour of the particle is investigated. For this 
purpose skewness 5* and kurtosis 7 are calculated. The non-Gaussian behaviour is clearly 
signalled in the nonzero value of these quantities. In general, the skewness is given by 

^ {{x-{x)f) ^ {x')-3{x){x') + 2{x)' 

and the kurtosis is calculated as 

_ {ix-{x)r) ^_ {x^)-A{x){x^) + 6{xr{x')-3{xr ^ 




FIG. 3: Mean square displacement of a spherical particle for (3RFq = 10 and different values of ( 



For the third and fourth moments of a: - in analogy to equation f|T4|) - one has to solve the 
integrals 

ft rt rt 



{{x{t)-xo)^) = 0^D^ / dti / &t2 / dig Fo3(cos0(ti) COS 0(^2) COS 0(^3)) 

Jo Jo Jo '- 

+ 3Fo(cOS0(ti))(/(t2)/(t3)) 



(22) 



and 

((x(t)-xor> 



Fq (cOS0(ti) COS 0(1(2) COS 0(1(3) COS 0(1(4)) 



/ dti / dt2 / dt3 / dt4 
Jo Jo Jo Jo 

+6F2(cos0(ti)cos0(t2))(/(t3)/(t4)) + {fih)f{h)fits)f{U)) 



(23) 



respectively. Before solving the time-integrals over the first summands, the time correlation 
functions 



(cos 



01 cos 02 COS03)j^ 



>i2>t3 = I cos(0o)e 



-Dr(tl-t2+t3) 



+ L-Driti+3t,-4t,) [cos(0o)e-^'*=' + cos(30o)e-9^^*^] (24) 
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and 



(cos 01 COS 02 COS 03 COS (j^A.) t,>t2>t3>U 

= ie-^^(*^-*^+*^-*^) [l + cos(20o)e-^^^**] 

^l^_^,(i,+3t,-4t3)| lg-D,(t3-t4) [1 + cos(20o)e-^^^*^] 

^ 1 -9D.(t3-t4) [cos(40o)e-^^^^*^ + cos(20o)e-^^'**] | (25) 

have to be evaluated. Here the notation 0j = 0(tj) with i G {1,2,3,4} is used again. Both 
in equation (l22l) and in equation ( |23l) . the remaining terms can be calculated easily using 
the expressions already obtained for the first and second moments. The complete analytical 
results for the third and fourth moments (and for the skewness and kurtosis) are presented 
in the appendix. 
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FIG. 4: Skewness S{t) for I3RFq = 10 and different values of 0o- 

Figures H] and [5] display the skewness 5* of the probability distribution of the particle 
position for different values of the initial angle 0o and the dimensionless quantity (3RFq 
which determines whether the self-propulsion or the motion due to the interaction with the 
solvent molecules is dominant. Figure H] shows that the sign of the skewness depends on 0o. 
If the x-component of the initial orientation is positive (— O.Stt < 0o < O.Stt), the skewness is 
negative, while initial angles between O.Stt and I.Stt lead to positive S. For symmetry reasons 
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FIG. 5: Skewness S{t) for (j)Q = and different values of I3RFq. 

the skewness is zero for (pQ = O.Stt. Further analysis of formula (H2i) (see the appendix) gives 
the leading long-time behaviour of the skewness S{t) as 



cos (0o) ((cos (0o))^ - 3) 



(26) 



3 (3 + 2a2) V4a2 + 6 
where the abbreviation a = [3RFq is used, i.e., the skewness decreases proportionally to 
t"^/2. Similar analysis of formula (jUj) for the kurtosis 7(t) reveals a long time behaviour 

as 



lit) 



-21a4 



(At) 



-1 



(27) 



9 -M2a2 + 4a4 

First of all, as can be seen from this formula and in figures [6HH1 the kurtosis does not depend 
on 00 for long times. The long-time tail, being proportional to is more pronounced than 
that for the skewness. Moreover, as displayed in figures [6] and [8|, for initial angles 00 7^ O.Svr 
the distribution is leptokurtic (positive kurtosis) for relatively short times and platykurtic 
(negative kurtosis) for relatively long times. Thus for intermediate times a change of sign 
is induced such that the kurtosis approaches its asymptotic value zero from below. This 
is in contrast to passive ellipsoidal particles in two dimensions where non-Gaussian 
behaviour is due to dissipatively coupled translational and rotational motion. In the latter 
case, the same scaling of the long-time tail proportional to l/t is found for the kurtosis but 
it approaches zero from above. 
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FIG. 6: Kurtosis ^{t) for (3RFq = 10 and different values of 
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FIG. 7: Kurtosis 7(t) for (pQ = and different values of (3RFq 



We expect that the different sign is linked to the one-dimensionality of our model rather 
than to the qualitatively different translation-rotation coupling, which is due to the driving 
force in our model as opposed to the different transverse and parallel short-time translational 
diffusivities in the passive ellipsoidal particle model. In particular, we expect the negative 
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FIG. 8: Kurtosis ^{t) for (pQ = O.Svr and different values of PRFq. 



kurtosis at long times t ^ D^^ to reflect a broad translational van Hove function 33|] with 
shorter tails as compared to a Gaussian distribution, which is attributed to the non-linear 
cos-term in equation ([1]). 



IV. RESULTS FOR A TIME-DEPENDENT TORQUE 

Let us now assume an additional internal or external torque. Before considering the case 
of an arbitrarily time-dependent torque M(t), we first consider a constant torque M. Solving 
the Langevin equations ([!]) and (jlj) under this assumption, one obtains 

{(f){t)) = 4>o + l3D,Mt = ^o + ujt (28) 

with the frequency u = fSD^-M and 

- m))f) = (29) 

By replacing 0o in formula ([9]) with 0o + ^t, the updated probability distribution of is 
gained. The mean position is obtained as 

{x{t) - xo) = Fo [D, cos(0o) - sin(0o) 

+e''^'\uj sin(0o + ^t) - D, cos(0o + ^t))] . (30) 

12 



In figure [9] this result is plotted for different values of the dimensionless quantity which 
is the ratio of the external torque over the thermal energy. The long-time mean position is 
given by 



lim {x{t) — xq) 



-Fq [A cos(0o) -uJsm{(t>o)] 



while the behaviour for short times is the same as in formula ( |T2l) for a vanishing torque. 
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FIG. 9: Mean position of a spherical particle with additional constant external torque for PRFq 
10, (/)o = and different values of (3M. 



Following the notation introduced in formula (iTSll the Green function is now given by 

1 / (01-02-^(tl-t2))' 



^AnD,{ti - t2) 



exp 



(32) 



This leads to 



(cos 



0lCOS02)t^ 



>t2 



-Dr(t,-t,) [cos(^(^^ - t2)) + COS (200 + toih + ts)) e-^^^*^] (33) 
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and by integration one obtains 

{{x{t)-Xor)=2Dt + (3'F,'D'i^ — 



D,t 



-DA 



2 + CU2 (D2 + ^2)2 



[ct;2 — cos{ujt) + 2ujD^ sin(co't) 



e~^'* ((-3D2 + tu^) cos(20o + ojt) 



(9Z^2+CJ2)(D2+CU2) 

+4L'i.t^ sin(20o + ^)) - + ^'^) cos(20o) - 4At^ sin(20o) 

1 

-lOL'rt^sin(20o + 2ojit)) - (12L)2 _ 2t^2^) cos(20o) + lOA^sin(20o) k34) 



e-^^^*((12L),2 - 2cu^) cos(20o + 2cjt) 



The result is displayed in figure dUl In this case, the long-time diffusion coefficient is given 
by 



D, = D + 



8 {PFpR^fP, 

9 (1 + (/5M)2) 



(35) 




Drt 

FIG. 10: Mean square displacement of a spherical particle with additional constant external torque 
for (3RFo = 10, (/>o = and different values of (3M. 



To generalize the preceding considerations, the torque M{t) is assumed to be arbitrarily 
time-dependent now. Similarly to the two special cases investigated so far, it can be seen 
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that the mean position of the particle is given by 



{x{t)) = (3FoD / cos 



'0 + /3A / M{t2)dt2 







e-^'-*Mti + xo. (36) 



The calculation of the mean square displacement starts with formula (fT4|) again. The first 
summand is the most interesting one because the other ones can be treated as before. Based 
on the formula 



(cos 



0icos02)i,>t, = ^e-^'(*^-*^)[cos(^/3A^'M(t)dt^ 

+ cos (^2(f)o + 2/3D, ' M{t)dt + /3 A ^ ' M{t)dt^ e'^^''^ (37) 



we introduce 







ut^ -.= (30, r M{t)dt. (38) 
Jo 

Using this notation the problem can be solved in a similar way as for a constant M. The 
mean square displacement is now given by 

{{x(t)-xof) = 2Dt + f3^F^D^ fdh r dt2e-^^(*i-*^) 

Jo Jo 

X [cositot, - ujt,) + COS (200 + + Wfe) e-^^^*2] . (39) 

V. CONCLUSIONS 



In conclusion, motivated by recent experiments on catalytic colloidal particles |15l . l23l . |24| , 
we have proposed and solved a model for a self-propelled colloidal particle on a substrate. An 
internal or external time-dependent torque is also included in the most general version of the 
model which can arise, e.g., from an external magnetic field. The first four moments of the 
particle displacement distribution were calculated analytically. Significant non-Gaussian 
behaviour was found for intermediate time. The normalized kurtosis changes sign and 
approaches zero from below with a massive long-time tail inversely proportional to time. 

Future work should address several generalizations of the model. First of all, the one- 
dimensionality of our model can be generalized towards higher dimensions both for the 
translational and orientational degrees of freedom. In particular, the translational degrees 

15 



of freedom can be considered to be two-dimensional (in a plane) , and the orientational ones 
on a sphere. For the latter case, first analytical results have been obtained Also, 
e.g., for weak gravity, the third translational dimension perpendicular to the substrate is 



getting important, which results in unusual sedimentation effects [35j. Furthermore, the 
self-propelled particle can be confined in the lateral direction 2J| which leads to a finite 
mean square displacement. This effect should be incorporated into a model study as well. 



36| and 



First results have been obtained for a circle- swimmer in planar circular geometry 
for swimmers in cuspy environments leading to self-rotating objects js? ]. 

Last not least, the collective behaviour of many interacting self-propelled par ticles is 



expected to lead to novel characteristic nonequilibrium effects both without 38|, |39|, |40|, |41| 



and with confinement 



371, 



42| . As stated in the introduction, the Smoluchowski equation. 



suitably generalized to self-propelled particles [421], is an appropriate starting point here and 



29|, 3^, l3l(| is expected to 



the general hierarchy of Bogolyubov-Born-Green-Kirkwood-Yvon 
be a valuable tool in order to derive approximations in a systematic way. This fact after all 
clearly links the present paper to the 100th anniversary of N. N. Bogolyubov. 
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APPENDIX 



Using the notation a = /SRFq and a scaled time r = D^t, we summarize here the analytical 
results for the third and fourth moments as well as for the skewness S and kurtosis 7: 

^3 ) = y arcos(0o) (1-e ) 

64 45 , , , 1 , , , 17 , , , _^ 

^27 " \~8 ^ ^ 24 ^ ^'^^ ^ y ^ 

1 1 
-— cos (3 0o) e"^ + 3 cos (0o) r - — cos (0o) e"^^ 
io 24 

+^ cos (3 0o) e"'" + ^ cos (0o) re~^ " ^ ^^'^ (3 0o)) , (40) 
16 



{x{t)-Xo)'^\ 64 2 256 



/ 3 



t' + '^ a'r (^e-^ + r - 1 + ^ cos (2 6) (e""^ - 4 e"^ + 3) 



261 1 19 1 

^16" + 600 ~ y + 240 

+— cos (4 0o) + — e"^^ - 1^ - — cos (2 0o) e"^^ 
192 ^ ^ ' 48 3 450 ^ ^ ' 

1 3 1 

cos (4 0o) e""^ + o cos (2 0o) cos (2 0o) e"^"^ 



120 ' ' 2 ' " ' 30 

229 1 ^ \ 

+— cos (2 0o) e-^ - — cos (4 0o) e"' ^ + 3 r cos (2 0o) e"^ j , (41) 



+ y (^"" + T - 1 + ^ cos (2 0o) (e-^^ - 4 e-^ + 3) ^ 



-3/2 



32 760 44 32 

x{ — — cos (0o) - 0? cos (0o) e"^ " 7^;; «^ cos (3 0o) e""^ - cos (0o) e"**"" 
9 ol Zl oi 

32 4 32 352 

a'^ cos (30o) e~'^'^ a^e"^"^ cos (3 Sn) H a? cos (30o) H cos (0o) re""^ 

135 ^ ' 405 ^ ^ 81 ^ ' 27 ^ ^ 

448 64 32 32 

cos (0o) e"^"^ + — a^e"^"^ cos (3 0o) - 7;;; a^e"^"^ cos (3 0o) - 77 a^e"^"^ cos (0o) 



27 



27 



27 



9 



8 
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